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In this article, we investigate the competing Glauber-type and Kawasaki-type dynamics with 
small- world network (SWN) effect, in the framework of the Gaussian model. The Glauber- type 
single-spin transition mechanism with probability p simulates the contact of the system with a heat 
bath and the Kawasaki-type dynamics with probability 1 — p simulates an external energy flux. Two 
different types of SWN effect are studied, one with the total number of links increased and the other 
with it conserved. The competition of the dynamics leads to an interesting self-organization process 
that can be characterized by a phase diagram with two identifiable temperatures. By studying 
the modification of the phase diagrams, the SWN effect on the two dynamics is analyzed. For 
the Glauber-type dynamics, more important is the altered average coordination number while the 
Kawasaki-type dynamics is enhanced by the long range spin interaction and redistribution. 
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I. INTRODUCTION 

Many systems in nature or society can be well described by small world networks (SWN), which were first proposed 
by Watts and Strogatz in Ref. [1]. SWNs are those intermediate between a regular lattice and a random network. 
They can be realized by introducing a very small portion of long range links to a regular lattice. Such networks 
appear as small- worlds like random graphs, i.e., with a small average node-node distance that grows logarithmically 
with the network size. Meanwhile they also have usually large clustering coefficients [2-4]. Since the first model 
system was proposed [1], a large literature has appeared and the properties of various models and processes on SWNs 
have been extensively studied, including percolation, coloring, coupled oscillators, iterated games, diffusion processes, 
epidemic processes, and spin models (see Ref. [2-4] and references therein). For example, it has been recently found 
that some spin-lattice models belonging to different universality classes now show mean- field behavior on SWNs [4-9] . 
Furthermore, recently Zhu et al. introduced SWN effect to critical dynamics [10], and the investigation has been 
extended to the kinetic properties of the spin models. 

According to Glauber's and Kawasaki's theory [11,12], the time evolution of the order parameters can be described 
by a Markov process with a specific choice of the dynamic mechanism. Two mechanisms have been extensively 
studied: Glauber's single-spin flip mechanism [11] with the order parameter nonconserved and Kawasaki's spin-pair 
exchange mechanism [12] with the order parameter conserved. In the past decade, an interesting problem has been 
attracting much attention, i.e., the competing Glauber- type and Kawasaki-type dynamics. We briefly explain the 
physical ground: Glauber's mechanism with probability p is used to simulate the contact of the system with a heat 
i-^J ■ bath and favors a low energy state. Meanwhile Kawasaki's mechanism with probability 1 — p simulates an external 
energy flux and, naturally, favors a high energy state. Then there will a competition between the two dynamics. 
O ■ This competing mechanism has been applied to the spin models [13-19], and interesting self-organization has been 
reported. All these works have considered regular lattice, and to extend the investigation to the small-world networks 
J> . should be of interest. At the same time, such a study may also have theoretical meaning. In Zhu's work [10], it 
has been found that the SWN effect may have different types of influence on the Glauber-type dynamics and the 
Kawasaki-type dynamics. By focusing on the competition of the dynamics, it is easy for us to further understand the 
■ influence of the SWN effect and highlight the disparities between the dynamics. 

In this article we investigate the SWN effect on the competing dynamics, and we demonstrate the conclusions by 
comparing the phase diagrams obtained with and without the SWN effect. Before we present the calculations, we 
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first describe the prototypes of SWNs that are used in the present study, and then briefly review the generalized 
competing mechanism. 

The SWN models: After the first prototype of SWN was proposed by Watts and Strogatz [1], there appear a 
variety of model systems that can be categorized into the two groups, i.e., with the total number of links increased 
or conserved. Correspondingly, in the present research we study two specific, while representative models described 
below: (1) in a one-dimensional loop, for example, each randomly selected pair of vertices are additionally connected 
with probability pa; (2) the vertices are visited one after another, and each of its links in the clockwise sense is left 
in place with probability 1—pR and is reconnected to a randomly selected other node with probability pr. Networks 
of higher dimensions can be similarly built. We call the first model adding- type small- world network (A-SWN) and 
the second one rewiring-type network (R-SWN). 

The generalized competing mechanism: Glauber's flip mechanism and Kawasaki's exchange mechanism were origi- 
nally presented for the Ising model. Recently they have been generalized to single-spin transition mechanism [20] and 
spin-pair redistribution mechanism [21], respectively, which can be applied to arbitrary spin systems. Correspondingly 
there is the generalized competing dynamics [19], which provides a basis for the discussions in this article. With the 
competing mechanism, the master equation can be written as 



where pG me denotes the single-spin transition with probability p and (1 — p) K me denotes the spin-pair redistribution 
with probability 1 — p. For your convenience we list in Sec. II and III some important equations of the transition 
mechanism and the redistribution mechanism, the details of which can be found in respective references. 

In this article we study the generalized competing dynamics in the framework of the kinetic Gaussian model, which 
is mathematically tractable. It has been revealed in Ref. [10] that the SWN effect on the Gaussian model is of 
the mean-field nature. Considering this, in the present study a simplified method is adopted, although a rigorous 
treatment is possible. 

In Sec. II we study the competing mechanics on A-SWN and In Sec. Ill, we extend the discussions to R-SWN. 
Section IV is the summarization with some discussions. 



In this section we present our discussion of the competing dynamics on adding-type small-world networks (the defi- 
nition see the Introduction). In subsection II A through II C we give the formulation of the Glauber- type mechanism, 
the Kawasaki- type mechanism and the competing mechanism, and then in subsection II D the competing mechanism 
is applied to the three-dimensional Gaussian model. 



With the Glauber-type single-spin transition mechanism [10,20] on A-SWN, the master equation can be written as 



^-P({a},t)=pG me + (l-p)K, 



(1) 



II. COMPETING MECHANISM ON A-SWN 



A. The Glauber-type mechanism 



£p(w,*) = -X;x;rwi(. 



cr 4 -> 6- t ) P ({a}, t) - W t (&i -> o-i) P ({(J^i), a, t)} , 



(2) 



i a 



where the transition probability 



Wi {<7i ^Vi) = Tr. ex P iWjjti] ' *»)] ' Qi = CXP \-P U ({ a j¥'} ' ^i)] ' 



With equation (2), one can continue to write that [10,20], 





where 




(4) 
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B. The Kawasaki- type mechanism 



With the Kawasaki-type spin-pair redistribution mechanism on A-SWN, the master equation can be written as 
[10,21] 



(jl) 
1 

-2 1 



+ \p a E E E [-^(^ w^w. *) 

+W,-|(CTj<7| -> CT J cr;)P({a J y i ,cr; 7 t/ s },<Tj,<T;,i)] , 



where the redistribution mechanism 



Wji(<7j<Ji — > (Tj-ct;) = — S^+a^aj+a, exp[-PH({a m } m7 tj,i, o-j, 07)], 



(5) 



With the master equation above, we can further get that 



q k (t) = Q% = -2Dq k (t)+J2J2 

{a} u 

+ PA {-(N -l)q k (t)+J2 



E VkWk,k+u (&k&k+oj — » 0"fc<7fc + u 



p(M,t) 



where denotes the summation taken over nearest neighbors. 



(6) 



C. Competing mechanism 

Naturally, both Glauber's dynamics and Kawasaki's dynamics favor a lower energy state. However, when the system 
under study is in contact with a heat bath while exposed to an external energy flux, one requires a competition between 
a process that favors lower system energy and another process that favors higher system energy. Usually, Glauber's 
mechanism is used to simulate the contact of the system with a heat bath and favors a lower energy state. Meanwhile 
Kawasaki's mechanism can be modified in order to simulate an external energy flux that drives the system towards 
higher energy. This can be achieved by switching j3 to — /?, or K = /3J = J/KbT to — K, and modifying the 
redistribution probability accordingly, i.e. 



Wji cx exp[-/3H] W jt oc exp[+j3H] 



(7) 



This means that the competition between Glauber's mechanism and Kawasaki's mechanism is actually that between 
ferromagnetism and antiferromagnetism. 

Based on the above considerations, we use the competing mechanism to simulate the dynamics of a system in 
contact with a heat bath and exposed to an external energy flux simultaneously. The master equation can be written 
as (Eq. (1)), 

d 



— P ({a}, t) = pGme + (1 - p) K me , 



where 



G me = -^^[W;(a 4 -» &i)P({a},t) - W t {&i - (Ti)P{{<T j¥ ^,a,t). 
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and 

+Wj;((Tj(Ti -> cr J -C7;)P({(7j 7 ti,cr; 7 t/ s },aj i (Ti,i)] . 
From the master equation, we can obtain the evolving equation of single spins, 

j t q k (t)=pQ<Z + (l-p)QZ. (8) 

The first term describes the Glauber-type dynamics which is used to simulate the influence of the energy flux and Q k 
is given by Eq. (3). The second term describes the Kawasaki- type dynamics which is used to simulate the influence 
of the energy flux. As mentioned above, we require it to favor higher energy state. However, the redistribution 
probability, Eq. (5), clearly favors a lower energy state. Considering this, we can use the expression of Eq. (6) for 
Q% , if we simply switch the sign before /? in the redistribution probability. 



D. The competing dynamics in the Gaussian model 

Now we turn to the Gaussian model built on an A-SWN, of which the Hamiltonian can be written as 

—ftH {{<?}) = if yVjjfc ((Ti+1JM + Ci,j+l,k + (Ti,j,k+l) 



ijk 
1 

2 1 



^ I 

+ ?PAKj2 (T i3k E Ui 'i' k ' ~ 2 PAK J2 a ijk' ( 9 ) 

ijk i'j'k' ijk 

where (3 = 1/fcsT and K = (3 J. We only discuss the case of K > 0, namely, J > 0. This case corresponds to the 
ferromagnetic system. The expression (9) is obtained by taking into consideration the fact that the influence of the 
system as a whole on individual spins is of the mean-field nature [10] (for related discussions see Ref. [4-9]). The 
spins can take any real value from — co to +oo. The probability of finding a given spin between a k and <r k + da k is 

assumed to be the Gaussian-type distribution, f((T k ) da k — \J1^ cxp (— |cr^) d<r k , where & is a distribution constant 

independent of temperature. Therefore, the summation of the spin value turns into an integration, and we can further 
obtain from Eqs. (3) and (6), 

K K 
Qfjk = -««*(*) + tE + U,j+u,k(t) + <li,j,k+u,(t)) + T (N - l)p A M(t), (K > 0) 



6 w =±l 



where M (t) = ±J2k Ik (i), and 



®^ jk ~ 2 [b+ (_jf)] k{[fa+M,fc _ H,j,k) - (li,j,k - Qi-i,j,k)] 

+ [(.Qi,j+l,k - Qi,j,k) - {li,j,k - Qi,j-l,k)] + [(Qi,j,k+1 - Qi,j,k) - {li,j,k - Qi,j,k-l)]} 

(-K) 

+ 2[b+ (-K)] ^ , ^ qi - 1 ' j ' k ~ 9i-i,j+i,k - Qi-i,j-i,k) + {2qi-i,j,k - Qi,j,k - 1i-2,j,k) 

+2{2qi+l,j,k - Qi+l,j+l,k - Qi+\,j-l,k) + (2<?i+l,j,fc - Qi,j,k - Qi+2,j,k) 
+2(2g ii j_i ) fe - qij-i, k+ i - qij-i t h-l) + {^Qi,j-i,k — Qi,j,k - Qi,j-2,k) 
+2(2q , j 1 j+l 1 fe — q-Lj + lM+l — Qi,j+l,k-l) + {^Qi,j+l,k — Qi,j,k — Qi,j+2,k) 
+ 2(2(?i,j,fc-l — <li+l,j,k-l — Qi,-lj,k-l) + {^Qi,j,-lk — Qi,j,k — Qi,j-2,k) 

+2(2g ijjfc+ i - q i+ ij tk +i - Qi-i,j,k+i) + {%Qi,j,k+i — Qi,j,k - Qi,j,k+2)] 
-M (JV - 1) [q ijk (t) M(t)] + PA (N 1)^1 

x [(»-H-,j,fc - M (t)) + {Qi,i+u,k - M(t)) + {qi,j, k+UJ - M{t)} , (b > K > 0) (10) 

o>=±l 
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where 6 > K > is required by the convergence of the integration. 

Now we determine the system behavior by studying the following order parameters. First, from Eq. (8) we obtain 



d_ 
~dt 



^E jM 1 ) - ^ E q& + a - f )^ E o«* 



6 



K 



-{N-1)pa 



M(t), 



(11) 



J2ijk Qijk = means that the Kawasaki-type dynamics does not change the value of M (t). Second, we define 



ijfc 



and similarly we obtain 



M 'w s ir E iafow = E o& + a - ^ E 

r/fc zjfc ijfe 



where 



and 



Q[% = -</ m (t) -f E + 3jj-H-,fc(*) + 9lj,fc-H-(*)] + f - l)PA(-l) i+J ' +fc M(t), 



(12) 



^ijk 



:&{[( l'i+l,j,k l'i,j,k) (l'i,j,k + 9i-l,j,fe)] 



2(6 -X) 

+ [( _ 9i,j+l,fe _ 9i,j,fe) _ (9i,j,fe + 9i,j-l,fe)] + [( _ 9i,j,fe+l _ 9i,j,fe) - fe'j.fc + 9i,j,fc-l)] } 
_ 2(6 — if) ^' i - 1 '3' k ~ + _ + (~^Qi-l,j,k ~ 5i,j,fe — 1i-2,j,k) 

+2(-2q- +lj fe - g^+ij+i^ - Q'i+i.j—i.fe) + (— 2<l'i+i,j,k ~ 9i,j,fc — 9i+2,j,fe) 

+ 2 ( _2 9i,j-l,fe _ 9i,j-l,fe+l - 5i,j-l,fe-l) + ( _2 5i,j-l,fe - Qi,j,k ~ 9i,j-2,fe) 
+ 2 ( — 2 <Zi,j+i,fc _ <7jj+i,fe+i _ <7jj+i,fe-i) + (~2g- j + i t k — ~ 9i,j+2,fe) 
+ 2 (— 2<li,j,k-l ~ 9i+l,j,fe-l _ <7i,-l?,fc-l) + ( 2qr ijJ - j _ lfc - (Zj^fc - Qij-2,k) 
+ 2 ( — 2 9i,j,fe+l - 9i+l,j,fc+l - 9i-l,j,fe+l) + ( — 2 5i,j,fe+l - 9i,j,fe - 9i,j,fc+2)] 

_^v-^{^ w _ ( _ 1)i+j+feMW 



Performing the summation over the indices i, j, and k, we get 

,6if-6 PA(iV-l) 



6-if 



6if 



M'(t). 



The Solutions of equations (11) and (13) are 

M(i) = M(0)exp j-p 

and 



(N - l) PA 



K 
K9 



M'(t)=M' (O)exp 



P 1 + 



if / 

i/ p = l,then K > 
i/f>^ l.tften (6 > if > 0) 



(1-P) 



36 



t|,(if >0) 

pa(JV-i) 



(13) 



(14) 



6 - K/K° c 



t 



(15) 



•5 



where K® = \ J\ /ksTc = 6/6 is the critical point of the three-dimensional Gaussian model without the SWN effect. 
(1) When the SWN effect does not exist and pa = 0, 



M(t) = M (0) exp [-p (1 - K/K° c ) t] 



(16) 



M'(t) = M' (0)exp 



P (l + K/K° c )+36(l-p) 



1 - K/K 



01 



6 - K/K° 



Analyzing the long-time system behavior, we find that: 
(l.a) For the case of 



K < K° c , (T > T c °), 

we have both vanishing M (t) and M' (t), which correspond to the paramagnetic phase, 
(l.b) For the case of 



M' (t) 0, and 



M(t) = M (0) exp 



K^K° 



t 

-P- 

T 



1 - K/K° 



oo. 



The critical slowing down of the order parameter M(t) will appear at the critical point K®. 
(l.c) For the case of 



K > Kl (T < T c ), 



and 



^ < kp (~ 36 + Alp + v^ 1561 ^ 2 - 2808 P + 1296 )) . 

we have nonvanishing M (t) and vanishing M 1 (t), which correspond to the ferromagnetic phase, 
(l.d) For the case of 



p (-36 + Alp + v/(1561p 2 - 2808p+ 1296)) < — ^ < 6 



(17) 



(18) 



(19) 



2p 

we have both nonvanishing M (t) and M' (t) . 

With both of the order parameters nonvanishing, this phase cannot be simply identified as ferromagnetic or anti- 
ferromagnetic. We name it as heterophase, and from Eqs. (16) and (17) we can see that the system behavior strongly 
depends on the initial condition and the temperature. 

(l.e) For the case of 

K/K° c > 6, 

we could obtain M (t) ^ 0, and when p = 1 (with pure Glauber- type dynamics simulating the contact with a heat 
bath), we have M' (t) — > 0. However, with p / 1 (with an external energy flow), because the antifcrromagnetic 
Kawasaki-type dynamics is limited by the condition (b > K > 0) to ensure the convergence of the integration, we 
cannot obtain the value of M' (t). As a result, we cannot theoretically obtain the system behavior in the temperature 
region K/K° > 6 1 . 

The phase diagram is shown in Fig. 1(a). 

(2) Now the A-SWN effect is introduced and pa ^ 0. We suppose pa = 1/JV, and then we have 



The same problem exists in Ref. [19], i.e., the system behavior cannot be theoretically obtained for the region K > b. 



6 



M(t) = M (0) exp <^ -p 



1\ K 



6 J K° 



(20) 



and 



M'(t) = M' (0) cxp |- (p (1 + K/K c ) + (1 - p) 



36 



6-K/RO 2 



(21) 



Analyzing the long-time asymptotic behavior, we can similarly get: 
(2. a) For the case of 

K < K c \ pA=1/N = Q -Kl 

the condition leads to vanishing M (t) and vanishing M' (t), which correspond to the paramagnetic phase. 
(2.b) For the case of 

A ^ ^c\ PA= i/ N — 

The critical slowing down of the system will appear. 
(2.c) For the case of 



79 + ^712 V - 12862p+ 5929 K 6 
2(3p-l) > K° > 7 

we have nonvanishing M (t) and vanishing M' (t), which correspond to the ferromagnetic phase. 
(2.d) For the case of 



K 89p - 79 + v/712939 2 - 12862p+ 5929 

we have both nonvanishing M (t) and M' (i), which correspond to the heterophase. 
(2.e) For the case of 

K/K° > 6, 

For the same reason (see (l.e)), we cannot obtain theoretically the system behavior. 
The phase diagram is shown in Fig. 1(b). 



III. COMPETING MECHANISM ON R-SWN 



Now we investigate the competing dynamics on rewiring-type small-world networks (the definition see the Intro- 
duction) , and the sequence of the content is the same as that in the previous section. 

With the Glauber- type mechanism, the master equation has the same form as Eq. (2), and the single-spin evolving 
equation is also given by Eq. (3). 

With the Kawasaki-type mechanism, the master equation can be written as 

j t P(W},t) = (l-Pfl)£E [-Wjiiajat - ^)P({a},i) 

(jl) 

+Wjl(&j&l — > (Tj(Tl)P({(Tjjti,aijik},0-j,0-l,t)] 

E E E i-wdw - w)p(Wh t) 

+Wji(&j&i -> (Tj<ri)P{{<rjTti,<rijtk},<7j,<ri,t)] ■ 
The redistribution probability Wji is given by Eq. (5), and the single-spin evolving equation 
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4«fc(t) = Qk = (1 - ph) <! -2l>g fc (t) + E E 

1 {CT }u;=±l 



+ 



Dp 



R 



N - 1 



-(iv-i) gfe (t)+E 



p(W,t)l. 



P({a},t) 



E E ^kW k i(a k (Ti -> <5- fe o- ; ) 



With the competing mechanism, the evolution of % (£) is once again described by an equation having the same form 
as Eq. (8) 

j t q k (t)= P Q% + (l-p)Q%, 

which consists of two terms, one corresponding to the Glauber-type dynamics and the other corresponding to the 
Kawasaki-type dynamics. In the following we turn to the Gaussian model. 

For the three-dimensional Gaussian model built on a R-SWN, the Hamiltonian can be written as 

—0H = K(l - p R ) E a ijk (cr»+l,j,fc + &i,j+l,k + Ci,j,fe+l) 
ijk 

3 PRK^2<T ijk E a i'3'k> ~ J^PR K ^Z^jk- 

ijk i'j'k' ijk 



N 



(22) 



In the single-spin evolving equation with the competing mechanism, the term that corresponds to the Glauber-type 
dynamics is now given by 



K 



Qfjk = -Qijk{t) + y (1 - p R ) E (n+^i,k + Qi,j+u,k + Qi,j,k+u) + -j^-KM(t), (K > 0) 



and the term that corresponds to the Kawasaki-type dynamics can be written as (similarly, we switch (3 to -(3, and 
K to -K) 



(1 K - 

^ijk ~ 



1 



2[b-K{l-p R )] 



b{[{qi+i,j, k - q%,j,k) - (Qi,j,k - Qi-i,j,k)] 

— {qi,j,k — qi,j-i,k)\ + \{qi,j,k+i — qtj, 

[2 (2qi_ij ik — qi-ij+i.k — + (2<7i-i,j,fc — qi,j,k — qi-2,j,k) 



2[b-K(l-p R )} 

+ [(qij+i,k — q%,j,k) — {qi,j,k — qi,j-i,k)\ + [(q%j,k+i — qt,j,k) — {qi,j,k — qi,j,k-i)]} 

K(l-p R ) 



+2 (2q. 
+2 {2q 
+2 {2q. 
+2 (2q. 



+ l,j,k — qi+l,j + l,k — <?»+l,j-l,fe) + {2qi+l,j,k ~ Qi,j,k — qi+2,j,k) 
,j-l,k — 9i,j-l,fe+l — qi,j-l,k-l) + (2<7i,j-l,fc — qi,j,k — 9i,j-2,fe) 

,j+i,k - q%,j+i,k+i - Qi,j+i,k-i) + (2<?»,j+i,k - qi,j,k - <7i,j+2,fc) 
,j,fe-i _ 9»+i,j,fe-i _ 5i,-ij,fe-i) + (2g , ij,-ife - g»,j,fc - qi,j-2,k) 



+2 (2qij tk+ i - q i+ ij^ k+ i - qi-i^M+i) + (2gi,j,fe+i - - qij^k+2) 



-^PR^qnk{t)-M{t) + 

(b>K(l-p R ) >0) 



tf(l-p*) 



E (li+u,j,k + qi,j+w,k + qi,j,k+u>) ~ 6M (i) 



.a>=±l 



where & > K{1 — p R ) > is required by the convergency of the integration. 

From the single-spin evolving equation with the competing mechanism, we can obtain 



^) = ^E^w=^-i + -j^ 



(23) 



and 
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|^'w=4e«& + (i-p)^Eq& 

ij K ij k 



= -p 



+(l-p){6 



.6K(l-p R )-b 3 Pl 



1 6K n 
b-K(l-p R ) 2 [ b 

The solutions of Eqs. (23) and (24) are, respectively, 

M(t) = M (0) cxp [-p(l - K/K° c )t] , (A > 0) 

and 



M'(i). 



M'(t) = M' (0) cxp 
+(l-p) 



1 + (l-^)^o 



36 



+ 



A' 

c 

3pr 



6-K(l-p R )/K° ' 2 

if p=l, then K > 
i/ p ^ 1, tften (6 > K{\ - p R ) > 0) 



l-(l-Pfl) 



A" 



where A° = J/fc B T c ° = 6/6 is the critical point without the SWN effect. 

(1) When no SWN effect is considered and p R = 0, the phase diagram is given by Fig. f (a). 

(2) Now we introduce the R-SWN effect and set p R = 0.1, and then Eqs. (25) and (26) become 



and 



M'{t) = M' (0) exp 



M(t) = M (0) cxp [-p(l - K/K° c )t] 



p i 



ioAoy 



36 



+ 



R 1 20 



9A \ 



1 io^cV 



Analyzing the long-time asymptotic behavior, we find that: 
(2. a) For the case of 

A<A C °, (T>T C °), 

we have both vanishing M (t) and M' (t), and this corresponds to the paramagnetic phase. 
(2.b) For the case of 

A - A c °, 

we shall observe the critical slowing down of the system. 
(2.c) For the case of 



A > A c °, (T < T c °), 



and 



A 25230p - 22230 + 15(V26017p 2 - 46842p + 21609 
AO < 2(621p-81) ' 

the condition leads to nonvanishing M (t) and vanishing M' (t), which correspond to the ferromagnetic phase. 
(2.d) For the case of 



25230p - 22230 + 150 V26017p 2 - 46842p + 21609 A 20 
2(621p-81) < AO < Y' 

we have both nonvanishing M (t) and M' (t) , which correspond to the heterophase. 
(2.e) For the case of 

A/A° > 20/3, 

For the same reason (see Sec. II. (l.e)), we cannot give theoretically the behavior of the system in this region. 
The phase diagram is shown in Fig. 1(c). 



IV. SUMMARY 



Comparing Fig. 1(a), (b) and (c), we find both similarities and differences. The phase diagrams have similar 
structures while the boundaries may be shifted by the SWN effect. In the following we briefly explain our observations 
and discuss the nature of the competing mechanism and the SWN effect. 

First we notice that the phase diagrams are all separated into four regions, with two special temperatures: one is the 
critical temperature K c — J/ksT c , another is the limit temperature A" max = J/ksT m i n required by the convergence 
of the integration. When the SWN effect does not exist, K c = K® = J/ksT® = b/2D is the critical point of the 
Gaussian model, where b is the Gaussian distribution constant and D is the space dimension, and K max = b (since 
b > K — \J\/ksT > is the region where we can assure the convergence of the integration in the Kawasaki- 
type dynamics of the antifcrromagnetic system). When the A-SWN effect is introduced, K c = 6+ ^ N 6 ^i^ PA K^ and 

^max = b = 2DK° (in the example studied D = 3, p A ~ 1/iV, and K c /K° = 6/7). When the R-SWN effect is 
introduced, K c — and if max = ^ PR ^ — rr~^\ ( m tnc example studied D = 3, pn = 0.1, and K max /K® = 20/3). 
Now we describe the phase behavior: Above the critical temperature T c = J/ksKc, one can only observe a disordered 
state, namely the paramagnetic phase, because of the dominating heat noise. Below critical temperature T c , the 
system begins to show some kind of order, and the system behavior is determined by the competition between the 
Glauber-type mechanism which favors a ferromagnetic state, and the Kawasaki-type mechanism which favors an 
antifcrromagnetic state. The result of the competition is subject to two factors, the probability of each mechanism 
and the temperature. However, the Kawasaki-type dynamics describing an antifcrromagnetic Gaussian system is 
limited by the condition T > T min , which is necessary for the convergence of the integration. As a result, in the region 
T < Tmin we cannot theoretically obtain the system behavior. 

Now we discuss how and why the two special temperatures T c and T m ; n may be affected by the SWN effect. When the 
Gaussian model is built on a R-SWN, the temperature T c (= T c °) remains the same. This is because, as is clear in all 
of the phase diagrams, T c is determined by the competition between the heat noise and the Glauber-type mechanism, 
which favors a lower energy state. One important characteristic of the Gaussian model is that this temperature can 
be further determined by the average coordination number. For example, without the SWN effect this temperature 
is given by J/ksTc = b/2D, where 2D is the average coordination number of a D-dimensional lattice. On a R-SWN, 
the average coordination number is not changed when a portion of the regular links are rewired, and thus T c also 
remains unchanged. On the other hand the temperature T m ; n is lower than the value on a regular network. On 
R-SWN long range links have made long range spin-pair redistribution possible, and therefore the influence of the 
Kawasaki-type mechanism is enhanced. Thus it is not difficult to understand why on R-SWN there is a lower value 
of the temperature T m ; n . 

When the Gaussian model is built on an A-SWN, the temperature T c 6 +( w ~ 1 )p a yo ^ becomes higher. As is 

explained above, this temperature is determined by the average coordination number, which is | • 2D on a three 
dimensional A-SWN with pa ~ 1/N. At the same time the temperature T m i n is unchanged. We may explain that it is 
because the two dynamics are both strengthened, i.e., by the larger coordination number and the long range spin-pair 
redistribution, and the temperature that characterizes a certain counterbalance remains the same. 

To summarize, in this article we investigate the competing Glauber-type and Kawasaki-type dynamics on two typical 
three-dimensional small- world networks, adding- type (A-SWN) and rewiring -type (R-SWN), in the framework of the 
Gaussian model. We get the evolution of the order parameters, M{t) and M'(t), and by analyzing the long time 
asymptotic behavior we draw the phase diagrams. With the competing mechanism, there exist two easily identifiable 
special temperatures. The influence of the long range links is analyzed: For the Glauber-type dynamics, more 
important is the altered average coordination number while the Kawasaki-type dynamics is enhanced by the long 
range spin interaction and redistribution. 
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Caption of figures 



Fig. 1 The phase diagrams of the three-dimensional Gaussian model on (a) regular lattice, (b) A-SWN, and (c) 
R-SWN. The regions Para, Fcrro, and Hctcro correspond to the paramagnetic, ferromagnetic, and heterophase phase, 
respectively, while the part of oblique lines belong to unknown region. 
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